Nonequilibrium conductance of asymmetric nanodevices in the Kondo regime 
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The scaling properties of the conductance of a Kondo impurity connected to two leads that are 
in or out of equilibrium has been extensively studied, both experimentally and theoretically. From 
these studies, a consensus has emerged regarding the analytic expression of the scaling function. 
The question addressed in this brief report concerns the properties of the experimentally measurable 
coefficient a present in the term describing the leading dependence of the conductance on eV/Tk, 
where V is the source-drain voltage and Tk the Kondo temperature. We study the dependence of a 
on the ratio of the lead-dot couplings for the particle-hole symmetric Anderson model and find that 
this dependence disappears in the strong coupling Kondo regime in which the charge fluctuations 
of the impurity vanish. 
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Universality of scaling plays a central role in the Kondo 
effeclii which describes the interaction of a magnetic im- 
purity with conduction electrons. As the temperature T 
is lowered below a crossover scale, denoted by the Kondo 
temperature Tk, the impurity spin becomes screened by 
conduction electrons. Recent advances in nanofabrica- 
tion techniques now allow for the experimental explo- 
ration of Kondo physics by attaching two conducting 
leads, which we denote left L and right R, to a smaller 
region, such as a semiconductor quantum dot (QD) or 
a single molecule transistor (SMT), each of which may 
act as an effective impurity spin. One typically applies a 
voltage difference V between the leads and measures the 
differential conductance G(T, V) = dl/dV which grows 
at low T and V due to Kondo correlations. This be- 
havior is expected to be described by a scaling form 
G(T, V) = GoF{T/Tk, cV/Tk) where Go = G(0, 0) and 
F is a universal scaling function. The scale Tk can vary 
between lOOmK in QD devices^i^ and 150K in SMTs.''-^ 

Such devices are more accurately described by the An- 
derson model (AM) which takes into account charge fluc- 
tuations of the impurity. The Kondo model (KM) is re- 
covered from the particle-hole symmetric AM in the limit 
[/ — > cx) where U is the charging energy in the AM. For 
both models, for low energies T, eV <^ Tk the leading 
corrections to the conductance are given by (up to a re- 
definition of Tk ) 
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with various values predicted for the coefficient a from 
the ky\§rl and Kondo model KM^^i^ Recently, a was 
measured by two experiments, one done by Grobis et. 
al. on a quantum dot device^^ and the other done by 
Scott et. al. on an ensemble of single molecule transis- 
tors,-^ where the QD and SMTs were tuned to the Kondo 
regime T,eV <^ Tk- In the QD experiment, Tk varied 
from 150 to 300mK by varying the gate voltage in a sin- 
gle device, and a value aqp — O.I±0.015 was measured. 
In the SMTs experiment, Tk ranged from 34 to 155K in 
29 different devices, and a showed a systematic deviation 
from the QD value, asMT = 0.051 ± 0.01 (see Refs. Ulil 



for the precise fitting range). Various possible explana- 
tions were pointed out for the systematic difference of a 
in SMTs. Among those, the relative asymmetry of the L 
and R coupling [denoted by A in Eq. ^] was considered 
as a relevant issue 

With this experimental motivation, we calculate a for 
arbitrary device asymmetry. We consider the particle- 
hole symmetric Anderson model (SAM) which generi- 
cally includes charge fluctuations in the impurity. We 
find that a is independent of the degree of L-R asymme- 
try only in the Kondo limit. Once charge fluctuations are 
included, there is a dependence of a on the L-R asym- 
metry. We compare our result to previous theoretical 
results and also comment on the relevancy to the experi- 
ments. We find that the low value of a measured in SMTs 
can not be accounted for within the symmetric Anderson 
model. 

Our phenomenological approach consists of a modifi- 
cation of Nozieres Fermi liquid (FL) theoryiS to account 
for charge fluctuations in the SAM. We pay special at- 
tention to the effect of shifting the Kondo resonance at 
finite voltage [see Eq. ([15])]. Our result is a generaliza- 
tion of Oguri's calculation of the conductance which used 
a non-perturbative result for the Green's function of the 
SAMj^ and is found to reduce to Oguri's result for the 
special case he considered. 

The model of a single Anderson impurity connected to 
L and R leads is 



H — Ha + Hd + Ht, 

k(7 i—L,R 



Ht^Yl H {dUkat + H.c.) , 

ka i—L,R 



(2) 



(3) 



where da- annihilates an electron with spin a in the quan- 
tum dot d-level, Ckai annihilates a conduction electron 
with momentum k and spin a in the i — L, R lead. 
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ek = hvpk, and vp is the Fermi velocity. In the SAM 
that we will consider here, — —U/2. It is convenient 
to define the L-R asymmetry parameter 



L — R asymmetry 



A = 



(4) 



The chemical potentials /i^ (i = L,R), satisfying fi^ — 
fiR = eV, are measured relative the the Fermi level de- 
fined at zero voltage. Then the ratio 



(5) 



describes the relative voltage drop across the L and R 
tunnel junctions which could depend on the capacitive 
couplings of the leads and QD or SMT and which we 
treat as a second L-R asymmetry parameter. 

One can define the retarded d Green's function as 
G^(e) = -i dte^''\d{t)S + dU{t)). With these def- 
initions, the current is given by the Meir-Wingreen for- 
mulaic 



r de[fL - fu]^ [-ImG«(e) 

-/ — OO 



(6) 



with f, = /(e - ^^^) = L,R), /(e) - [1 + exp(e/r)]-i, 
A = + Ffl , Fi = T:vvf {i — L, R), and ly is the density 
of states in the leads. 

Below, we formulate an effective theory for the SAM 
that allows us to obtain the Green's function G^{e) for 
e, T, eV ^ A [A is the characteristic energy scale defined 
more explicitly below Eq. ([8]). In the Kondo limit A 
Tk\- We define linear combinations of the annihilation 
operators for the L and R leads 



(s) _ VLCkaL + VBCkaR 



(p) _ -VRCkcrL + VLCkaR 



(7) 

where v = \/\vl\^ + \vji\-^. Only the s— wave particles are 

coupled to the d-level, since Hf = J2ka- ^'^t'^'kl + H.c. 

The notion of a local Fermi liquid (FL), due to 
NoziereSjiS, was originally applied to the Kondo model 
but is actually more general and can be applied to the 
AM and, in particular, to the less complicated SAM. 
The quasiparticles of this FL theory are simply scatter- 
ing states whose incoming part coincides with that of the 
the s— wave particles a\f^ [the precise definition is given 
after Eq. (ITOl) ]. The theory itself consists of a low energy 
expansion of their scattering phase shift as a function of 
energy e (measured from the Fermi level) and of quasi- 
particle density Ua-, 



da ^ Oo + ^ ^ 

A j/A 



(8) 



where f —i, i —tj ^ is the energy scale over which the 
phase shift varies in the SAM, and /3 is a coefficient to 
be determined. For the general AM, the scattering phase 
shift at the Fermi energy So can be extracted using the 



Friedel sum rulei^i 6o = Im[ln G-^(0)|T=ey=o] — com- 
bined with exact results for G^. In the SAM, particle- 
hole symmetry and the adiabatic connection to the U = 
case implies Sq — tt/2. 

The Wilson ratio R — {5x/x)/{^Gv/Cv) is the ratio 
between the relative impurity contribution to the sus- 
ceptibility and to the specific heat. It can be calculated 
from the phase shift expansion, Eq. ([5]), to bei^ 



R=l + I3. 



(9) 



We will use this equation to determine 13 in terms of R 
which is a parameter describing the amount of charge 
fluctuations in the SAM. 

As an equivalent way to determine /3, consider an en- 
hancement of the Fermi energy by an amount e by adding 
to the Fermi sea a density of quasiparticles rio- = (j^e). 
For the KM, Nozieres argued that at energy e, corre- 
sponding to the new Fermi energy, one has 5a — 5^ since 
the Kondo resonance is tied to the Fermi level. Using 
Eq. ([8]) and Ua = vt^ this implies /? = 1. This argument 
should be modifled for the AM: a shift of the Fermi level 
by this transformation also implies that ed — >■ Cd — e, as 
measured relative to the new Fermi level at -|-e. There- 
fore a finite amount of charge ebud enters into the d-level 
which, for small e, is given in terms of the charge sus- 
ceptibility 8nd = — -^e- Here rid = Due to 
the Friedel sum rule,!* the phase shift at the new Fermi 
energy after this transformation is different than 6q and 
given by (5, = 7rnd/2 = - (7r/2)^e. Using Eq. (0) 

and ria = vt implies ^ = \ + 

This phase shift expansion can be equivalently de- 
scribed by a Hamiltonian 



SH 



Hoia<-P^) + Ho{b) + SH, 
27ri^A ~-~ 

kk a 



(10) 



where Ho{ip) = J^ka '^ki'la'^^ka (ip = a'-P\b). This Hamil- 
tonian describes the two last terms in the phase shift 
expansion Eq. 

The first term in Eq. ([5]), Sq, is incorporated into 
the definition of the 6-particles in Eq. (fTO|) . These fa- 
particles are single particle scattering states that de- 
scribe an incoming s-wave suffering a scattering phase 
shift Sq at the boundary. Formally, to define the fa- 
particles, one uses the unfolding transformation^^ where 



^^J'{x) = (l/27r) / d/ce-*'=^4"j {x G {-oo, cx^}) is a chiral 
field describing an s-wave scattering state with the left 
moving convention such that 2; > is the incoming part 
and a; < the outgoing part, x = being the bound- 
ary. From the definition of the phase shift Sq we have 
^i'\0~) = e2^*oV^-'^(0+). We define the fa-particles in 
terms of a scattering state with Sq — 7t/2 

^W(a;)=V^^)(x)sgn(x), (11) 
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and its Fourier modes bka = J dxe^'^^ ^'^\x) . 

Now consider eV ^ 0. As long as eV <^ A the system 
remains in the vicinity of the fixed point and the state 
at finite eV can be treated within the FL theory as a 
state with a non-thermal distribution of quasiparticles. 
We first consider single particle scattering states incom- 
ing from lead i = L,R. In second quantization those 
particles are annihilated by [ckaiY"' ■ The occupation 
of those incoming waves is simply {{c\^iy^{ck'iy'i')™) = 
Skk'Saa'Sw fi{ek)- Using Eq. ([7]), and the fact that before 
the scattering region (x > 0) the wave function of states 



and bka coincide, we have^^ 



(CkaL)"' = {vLbkcr ~ VRa-^})/v, 
{CkaR)"' - {vRhka + VLa^^})/v. 

This gives the nonequilibrium distribution function for 
the 6-particles 

{bla^k'a') = Skk'Saa' [TLfdek) + Ffl/fllefc)] /A. (12) 

Since the occupation of the fe-particles differs from 
the one defined at T — eV = 0, the second 
term of 6H generates a constant elastic scattering 
V h, where 



^te^'^* ■) = ^(ri/^L + THf,R)/A. (13) 



Fermi energy is given by [see Eq. 



(14) 



where fix is a nonequilibrium shift of the resonance, 



(15) 



This shift can be nonvanishing if eV ^ 0, and J7 > (or 
equivalently i? ^ 1). In the Kondo limit, = 2, Eq. (fTS)) 
implies that the resonance position [Ik shifts together 
with the average chemical potential \i — (/^^ -I- ^ji)/2: 
under /2 — -I- (5/i at fixed eV , Eq. ((T5|) gives ^.k 
UK + SfJ.. 

In order to calculate the current using the phenomeno- 
logical Hamiltonian Eq. ([TU|) . one can use the Meir- 
Wingreen formula, and relate the Green's function G^ (e) 
to the s-wave T-matrix, r^(e) = vG^{e)v. One can de- 
fine a T— matrix, T, for the 6-particles due to SH. It has 
an inelastic part denoted by Ti„. The relation between 
Ts and T, accounting for the small inelastic term T^, 
readsi^ 



1 

2i 



1 



„2i5s(e) 



-TTvfinie) 



As a result, the phase shift at energy e relative to the 



(16) 

The leading contribution to Tin originates from the 
diagram shown in Fig. [TJ containing three propagators 
of ^-particles whose occupation is given by Eq. (|12p . We 
find, using the Keldysh technique. 



where i(e) - 1 - 2 ^^/^(e) + r^/flXe)] / A and k = (ri/ii-Hr^^^)/[A(AiL -Mi?)] - (BF^ - rfl.)/[(l + B)A]. Plugging 
Eqs. and pT)) into Eq. and using the Meir-Wingreen formula Eq. ([B]), with G^{e) = v^^Ts{e), we obtain the 
conductance in the scaling form of Eq. ([TJ with Go = ct = [^+^(-^~^) 1 ^ Tk A and 

" = 2^ ' '^TT^B + 3 (i + j3)3 + - + j / (1 + 2(i? - 1)^) . (18) 

I 



The coefficient a can be expressed as a function of 3 
independent variables such as the Wilson ratio R and 
the L-R asymmetry parameters A and R. In general, a 
depends on the the L-R asymmetry parameters however, 
in the strong coupling limit C/ 00, equivalent to i? — > 
2, this dependence completely disappears from Eq. (jlSp 
and we obtain a\R = 2] = = 0.1519. In this limit, 
A — >■ 47r-\/u/27rexp[— 7r^u/8 + l/(2u)] which is the known 
expression for Tk, where u — [//(ttA). 



The value of a measured in QDs can be accounted for 
by charge fluctuations due to finite U since, for i? ^ 2, 
Eq. (Ull) gives a\R, A, B\ in the range 3/47r2 = 0.0759 < 
a < 0.3039 = 3/7r^. Since one can tune the gate voltage 
in QDs and this corresponds to tuning e^, it is plausible 
that the SAM applies at one value of the gate voltage cor- 
responding to maximal conductance. However the value 
asMT — 0.051 ± 0.01 measured in SMTs is lower than 
the expectation from the SAM. We note that a value 
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FIG. 1: Diagram used to calculate Ti„. Each line is a 
6-particle propagator, with the nonequilibrium distribution 
Eq. dUl). 



validity of our reported expression for a. 

As another application of Eq. pO)) . one can calculate 
the shot noise S = rft({(5/(t), «}) for the SAM to 

leading order in 1/A where SI = I — (I). At T — 
and A — B — 1, using results for S based on the ef- 
fective Hamiltonian SH of Eq. ([T0| with arbitrary coef- 

ficients,^^ we obtain S = ^ (^eV^j" ^ ^he 

ratio 



a = 0.157 ± 0.005 was also measured^a in Al/AlOx/Sc 
planar tunnel junctions at low temperatures^ and can 
be accounted for in our theory. 

We compare Eq. (|18p to the results of other ap- 
proaches for the KM and SAM. Firstly, our result is 
fully consistent with the results of Oguri^ based on 
Ward identities, however he concentrated on the spe- 
cial case oi A = B = 1; Kaminski, Nazarov, and Glaz- 
mani^ find Q!kng = gf^ for the KM for any A; Konik, 
Salcur, and Ludwigi find oksl = 4/7r^ for the SAM 
for j4 = 1 and large J7/A; Pustilnik and Glazman^- 
find for the KM for A < 1, apc = Rincon, Ali- 
gia, and Hallberg^ studied the SAM for the case B = 
^r/^l- a mistake was found in Eq. (10) in their pa- 
per, whereas the corrected formula^^ ^^In-La — 1 — 



ii) 

with our Eq. ([TSl 



7r''(l+2(J?. 
3 



4-3A+(2+3A)(J?.-l)' 



(f) 



agrees 



Given that aKSL, which differs 
from our result, was obtained approximately and is not 
claimed to be exact, and given the agreement with the 
exact formulation due to Ogurij^ we are convinced of the 



2h 



l-|-9(i?-l)2 
l + 5(i?- 1)2 



e, 



(19) 
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where = — I, can be interpreted as a backscatter- 
ing charge.— This charge crosses from e* = e for i? = 1 
(noninteracting resonance level) to e* — 5/3e for R = 2 
(Kondo resonance))^"— 

In conclusion, we extended Nozieres Fermi liquid the- 
oryiS to account for charge fluctuations in the parti- 
cle hole symmetric Anderson model and calculated the 
transport coefficient a present in the term describing the 
leading dependence on eV/Tx- We included explicitly 
the effects of L-R asymmetry of the device and discussed 
the relation to recent experimental results J^i^ 

After this work was essentially completed, we became 
aware of another work^^ that obtains a = for the 
KM with arbitrary L-R asymmetry, consistent with our 
result in the special case without charge fluctuations. 

We thank I. Affleck, L. I. Glazman, D. Goldhaber- 
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This work was supported by NSERC. 



^ J. Kondo, Prog. Theor. Phys 32, 37 (1964). 

^ D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, 
D. Abusch-Magder, U. Meirav, and M. A. Kastner, Na- 
ture 391, 156 (1998). 

^ S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwen- 
hoven, Science 281, 540 (1998). 

J. Park, A. N. Pasupathy, J. I. Goldsmith, C. Chang, 
Y. Yaish, J. R. Petta, M. R. amd J. P. Sethna, H. D. 
Abrua, P. L. McEuen, and D. C. Ralph, Nature 417, 722 
(2002). 

^ W. Liang, M. P. Shores, M. Bockrath, J. R. Long, and 
H. Park, Nature 417, 725 (2002). 

® T. A. Costi, A. C. Hewson, and V. Zlatic, J. Phys.: Con- 
dens. Matter 6, 2519 (1994). 

^ R. M. Konik, H. Saleur, and A. Ludwig, Phys. Rev. B 66, 
125304 (2002). 

* A. Oguri, Phys. Rev. B 64, 153305 (2001). 

® J. Rincon, A. A. Aligia, and K. Hallberg, Phys. Rev. B 79, 
121301 (2009). 

° A. Schiller and S. Hershfield, Phys. Rev. B 51, 12896 
(1995). 

^ K. Majumdar, A. Schiller, and S. Hershfield, Phys. Rev. B 

57, 2991 (1998). 
^ M. Pustilnik and L. I. Glazman, J. Phys.: Condens. Matter 

16, R513 (2004); ibid, Lectures notes of the Les Houches 



Summer School 2004, cond-mat/0501007. 

A. Kaminski, Y. V. Nazarov, and L. I. Glazman, Phys. 

Rev. B 62, 8154 (2000). 

M. Grobis, I. G. Rau, R. M. Potok, H. Shtrikman, 
and D. Goldhaber-Gordon, Phys. Rev. Lett. 100, 246601 
(2008). 

15 G. D. Scott, Z. K. Keane, J. W. Ciszek, J. M. Tour, and 

D. Natelson, Phys. Rev. B 79, 165413 (2009). 

1® P. Nozieres, J. Low Temp. Phys. 17, 31 (1974). 
" Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992). 

1® D. C. Langreth, Phys. Rev. 150, 516 (1966). 

^'^ See, for example, I. AfHeck and A. W. W. Ludwig, Nucl. 

Phys. B 360, 641 (1991). 

Juhn-Jong Lin, private communication. 
^1 S. S. Yeh and J. J. Lin, Phys. Rev. B 79, 012411 (2009). 

A. A. Aligia, private communication. 

A. O. Gogolin and A. Komnik, Phys. Rev. Lett. 97, 016602 
(2006). 

E. Sela, Y. Oreg, F. von Oppen, and J. Koch, Phys. Rev. 
Lett. 97, 086601 (2006). 

^5 T. Fujii, unpublished, arXiv:0807.4869vl. 

C. Mora, P. Vitushinsky, X. Leyronas, A. A. Clerk and K. 
Le Hur, Phys. Rev. B 80, 155322 (2009). 



